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Abstract. We study the statistics of quantum transmission through a one-dimensio- 
nal disordered system modelled by a sequence of independent scattering units. Each 
unit is characterized by its length and by its action, which is proportional to the 
logarithm of the transmission probability through this unit. Unit actions and lengths 
are independent random variables, with a common distribution that is either narrow 
or broad. This investigation is motivated by results on disordered systems with non- 
stationary random potentials whose fluctuations grow with distance. 

In the statistical ensemble at fixed total sample length four phases can be 
distinguished, according to the values of the indices characterizing the distribution 
of the unit actions and lengths. The sample action, which is proportional to the 
logarithm of the conductance across the sample, is found to obey a fluctuating scaling 
law, and therefore to be non-self-averaging, in three of the four phases. According to 
the values of the two above mentioned indices, the sample action may typically grow 
less rapidly than linearly with the sample length (under localization), more rapidly 
than linearly (superlocalization), or linearly but with non-trivial sample-to-sample 
fluctuations (fluctuating localization). 
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1. Introduction 

The phenomenon of localization of a quantum particle by a random potential is well 
understood (see, e.g., [1, 2]), especially in the one-dimensional case (see, e.g., [3]). 
It is known that this phenomenon has a major influence on coherent quantum 
transport. The amplitude T of transmission through a sample is therefore a quantity 
of central importance in this context. The interpretation of experiments about the 
effects of localization on transport is done with the help of the (two-probe) Landauer 
formula [4, 5], which relates the probability of transmission |T| 2 through a sample to 
its conductance g at zero temperature, 

In the one-dimensional disordered systems usually considered, namely those 
modelled by stationary random potentials with short-range correlations, all the 
eigenstates are exponentially localized. Thus any sample whose length L is much greater 
than the localization length £ is an insulator [1, 2, 3]. The probability of transmission is 
a strongly fluctuating quantity in the insulating regime and so it is preferable to consider 
its logarithm [6, 7]. It is therefore appropriate to introduce the sample action 

S = -\\n\T\\ (1.2) 

The statistics of the transmission probability in the insulating regime are determined 
by the distribution of S. The average of the sample action over all the configurations of 
the random potential grows asymptotically as 

(S) » |. (1.3) 

The sample action S is self-averaging in the strong sense that each cumulant of its 
distribution increases in direct proportion to the sample length [3]. Said otherwise, the 
distribution of S is similar to that of the free energy of a one-dimensional disordered 
thermodynamical system; the transfer matrix formalism [3, 8, 9] is well adapted to 
demonstrate this deep analogy. 

The situation is quite different for one-dimensional disordered systems modelled by 
non-stationary random potentials whose fluctuations grow with distance [10, 11, 12]. 
Typical eigenstates turn out to be superlocalized rather than localized in samples whose 
length exceeds a well-defined crossover length £ CT , which diverges in the weak-disorder 
regime. The action S for long enough samples (L 3> £ CT ) exhibits two novel features with 
respect to the usual one-dimensional localization problem: it grows faster than linearly 
with L and it is not self-averaging in the sense that it keeps on fluctuating in arbitrarily 
large samples. For instance, in the case of a self-affine Gaussian random potential whose 
fluctuations increase as a power law with a Hurst exponent H in the range < H < 1, 
the scaling law of the sample action is the following [12]: 

(1.4) 
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Here E denotes the energy of the particle, £ cr the above mentioned crossover length, 
and y a rescaled random variable whose distribution is a universal function in the sense 
that it is determined by one parameter only, namely the value of H . 

In this work we address the question whether the properties of the two previous 
classes of models can be combined into a single model of one-dimensional disordered 
system, in which the concomitant features of self-averaging and conventional localization 
hold in a certain range of parameters but are replaced by the concomitant properties 
of lack of self-averaging and departure from conventional localization in another range. 
More specifically, the purpose of this paper is to present and study a model which 
exhibits the above properties and is simple enough to be exactly solvable by analytical 
means. The article is organized as follows. In Section 2 we describe the model and 
set the notation used throughout the paper. In Section 3 we study the statistics of 
the transmission probability through samples consisting of a large but fixed number N 
of units. In Section 4 we study the statistics of the transmission probability through 
samples whose length L is large but fixed. The results and their physical consequences 
are discussed in Section 5. Finally, an appendix is devoted to the derivation of the 
formula of the sample action on which our calculations are based. 

2. The model 

The model considered in this paper consists of a disordered array of elementary units, put 
together end to end on a line. Each of these units consists of a disordered part connected 
to two perfect leads, as shown pictorially in Figure 1 where disordered segments are 
represented by boxes, leads by line segments, and the connection between each unit and 
the next one by a dot. 
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Figure 1. Schematic picture of the model considered in this paper, consisting 
of a disordered array of elementary units whose lengths l n and actions s n are 
independent random variables. 

The nth unit has a length £ n and its transport properties are characterized by the 
transmission and reflection amplitudes t n , r n , and r' n . We shall work in the regime in 
which the probability of transmission through every unit is very small {\t n \ 2 <^1). As 
explained in the appendix, in this weak transmission regime and with the simplifying 
assumption that the reflection phases are random and uniform, the law of addition for 
the actions (A. 24) holds. 

The basic observables of our model are therefore the following. For a sample made 
of N units, we consider its total action Sn and total length L N . These quantities are 
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simply expressed as sums, 

N N 

S N = Y,Sn, L N =J2 £ n, (2.1) 
n=l n=l 

where the actions s n and lengths £ n of the units are the microscopic variables of the 
model. The simplicity of these definitions of the basic observables makes the present 
model exactly solvable by analytical means. 

Disorder is introduced into the model by assuming that the unit actions s n and the 
unit lengths i n are mutually independent random variables. The actions s n are assumed 
to be greater than some minimal action s min > and to be distributed identically, with 
a common continuous probability distribution defined by the density f s (s). Similarly, 
the lengths £ n are assumed to be greater than some minimal length £ min > and to be 
distributed identically, with a common distribution defined by the density f e(£). Thus 
the sequences of random variables Sn and L N obtained by adding units to a sample 
constitute examples of a renewal process (see, e.g., [13]). 

The samples may be classified in two different ways, either by their number N of 
units, ignoring their length, or by their length L, ignoring the number of units they 
consist of. These two different points of view give rise to the following two different 
statistical ensembles. 

• Statistical ensemble of the samples with a fixed number N of units. In this ensemble 
the sample action Sn and the sample length Ln, equation (2.1), are the sums of a 
fixed number N of independent random variables. 

• Statistical ensemble of the samples with a fixed length L. In this ensemble the 
sample action, denoted by Sl, is the sum of a random number N L of independent 
random variables, 

N L 

Sl = Y, (2-2) 

71=1 

For a given configuration {£ n } of the unit lengths, Nl is defined as the number of 
units which are entirely contained in the sample of fixed length L. This choice is 
motivated by simplicity. The random number Nl therefore obeys the inequalities 

N L N L +1 

L Nl <L< L Nl +i, i.e., Y j i n <L< J2 £ n - (2.3) 

71=1 n=l 

Since the length of any unit is greater than the minimal length 4m n , the random 
integer Nl can take any value between and its maximal value [L/£ min ], the 
symbol [x] denoting the integer part of the number x. 

These two statistical ensembles will be carefully distinguished throughout the 
following and successively investigated in Sections 3 and 4. The statistics of the 
transmission may indeed exhibit qualitative differences between both ensembles. 
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Generalizing a notation already introduced, we shall denote the probability density 
of a random variable 77 by f v (v)- The moment generating function and moment function 
of i] will be denoted by f v (x) and m v (p) respectively. These functions are defined as 

Ux) = <e— ') = J e-* fM d77, m v (p) = (rf) = J if f v ( V ) drj, (2.4) 

where the real parts of the complex variables x and p take their values in ranges such 
that the integrals are convergent. The moment function is only defined for a positive 
random variable rj. The inverse formulae are 

fM = I S eXV ^ {x) = I^[ ^ l mr > {p) - (2 - 5) 

The moment generating function and the moment function are linked together by the 
identity 

I r°° ~ 

m ^-P) = Y(p)L ^ xP ~ l&x (R e P>0), (2.6) 

where the symbol Kez denotes the real part of the complex number z and T(z) the 
gamma function. The identity (2.6) is easily verified by using the definition of f v (x), 
equation (2.4), on the right-hand side and evaluating the integral in x with the help of 
the first of the pair of identities ([14, p. 317]) 

r(p) 



J™ x 1 *- 1 e~ xy dx = (Rep>0, Rey>0), 



dy e xy x 13 - 1 nA (2-7) 



f dy x' 

I — = — — — - (Rex>0). 

J 27ri v p r p) v ; 



27ri y p T(p) 

Rey>0 y KFJ 

It is convenient to characterize the probability densities f s (s) and fg(£) by the 
exponents /x, v and the scales sq, £0 of the power law fall off of the corresponding 
complementary distribution functions: 

/ s (si)dsi « ( — J (s -> 00), 



Prob > ^ = jf /^) d4 » ( j - 00). 

The notation x ~ y used here means that the ratio x/y tends to unity in the considered 
limit; it will be used extensively in the following. We shall also use the notation x ~ y, 
which has the weaker meaning that x and y are only proportional in the considered 
limit, up to an unimportant prefactor. 

We shall consider the following classes of probability distributions f s (s) and fg{£). 

• Narrow distributions (fi > 2; v > 2). These are the distributions whose first two 
moments at least are finite. Thus we have 

f s (x) = l-x(s)+ 1 1 x 2 (s 2 } + ---, 

f t (x) = l-x{Q + ±x 2 (P) + --., { } 

where (s) and (s 2 ) (resp. (£) and (£ 2 )) are the first two moments of f s (s) (resp. ft{£))- 
Said otherwise, the narrow distributions are those for which the central limit 
theorem holds. 
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The narrow distributions with a complementary distribution function that falls off 
more rapidly than any power law (we have then formally \x = oo or v = oo) have 
the property that all their moments are finite, and so 

fc>o \k>i K - j 

fe>0 /i - \fc>l K - 

where (s k ) and ((s k )) (resp. (£ k ) and are the moment and cumulant of order k 

oif s (s) (resp. f e (£)). 

Broad distributions (/j, < 2; v < 2). These are the distributions whose second 
moment is infinite. In other words, the broad distributions are those for which the 
central limit theorem does not hold. Two cases have to be distinguished. If ji < 1 
(resp. v < 1), the first moment of f s (s) (resp. fe(£)) is also infinite. On the other 
hand, if 1 < /i < 2 (resp. 1 < v < 2), the first moment of f s (s) (resp. fe(£)) is still 
finite. The expansions of the moment generating functions as x — > + begin as 

'i-T(i-fi)(s xy + --- (/i<i), 

v 1 - (s)x + |r(l - fi)\(s xY + ■■■ (l<n<2), 



fs(x) 



ft{x) 



(2.11) 

l-T(l-u)(£ x) v + --- (u<l), 



{ 1 - (£)x + |r(l - v)\%xf + ■■■ (1< v < 2), 
where the dots represent regular and singular terms of higher order [15]. Disorder 
modelled by broad distributions will be called broad disorder. 



3. Transmission statistics in the ensemble of samples with a fixed number 
iV of units 

The statistics of the transmission probability through samples with a fixed number iV 
of units are determined by the distribution of the sample action Sn, equation (2.1). 
We shall be interested here in the scaling properties of Sn in the large iV limit. This 
study may be considered as an introduction to that done in Section 4, in which the most 
interesting results are derived. 

It is useful to introduce the moment generating function fs N (x) of Sn- Since the 
unit actions are mutually independent random variables with a common distribution, 
fs N (x) takes the form 

fsA-r) Ur) X - (3.1) 

The scaling law of Sn can be obtained from the expression of this quantity in the 
large N and small x limits. Three cases have to be distinguished, as the expansions (2.9) 
and (2.11) of f s (x) for small x are different for /i > 2, /i < 1, and 1 < /i < 2. We shall 
discuss these three cases below in this order. The results of this section also apply 
to the distribution of the sample length Ln, provided we replace the index // and the 
parameters (s) and s in the obtained results by the index v and the parameters (£) 
and £ respectively. 
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3.1. Self-averaging with normal fluctuations: // > 2 

This is the easiest case to study, as the distribution of the unit actions s n is narrow. 
Using the expression (2.9) of f s (x) in (3.1), we obtain 

f SN (x) » exp (N (-x(s) + ±:r 2 (( S 2 »)) (3.2) 
in the relevant regime (N large and x small). Setting 

S N ^N(s) + (N((s 2 ))) 1 / 2 r 1 , (3.3) 
we find that the moment generating function of the rescaled random variable i] is 

Uv) = e y2/2 - (3-4) 

Equation (3.3) gives the scaling law of the sample action in the present case. The 
fluctuations of Sn about its mean N(s) increase as iV 1 / 2 , and so become negligible in 
relative value in the large N limit. Hence Sn is self- averaging. The distribution of rj is 
Gaussian, 

UV) = ■ (3.5) 



In the case of narrow distributions with a complementary distribution function that 
falls off more rapidly than any power law, S n is self-averaging in the strong sense that 
each cumulant of its distribution is exactly proportional to N: 

((S N )) = N((s k )). (3.6) 
This can be checked by using the expression (2.10) of f s (x) in (3.1). 

3.2. Fluctuating scaling law: ji < 1 

In this case the distribution of the unit actions s n is broad, with an infinite mean (s). 
Using the first expression (2.11) of f s (x) in (3.1), we obtain 

fs N (x) « exp (-AT(1 - fi)(s xr) (3.7) 

in the relevant regime (N large and x small). Setting 

SN^Soiril-riY^N 1 ^^, (3.8) 

we find that the positive rescaled random variable £ M has a non-trivial limiting 
distribution, whose moment generating function is 

/c,(y) = e"^. (3.9) 

Equation (3.8) gives the scaling law of the sample action in the present case. It 
shows that Sn increases as N 1 ^, that is, faster than the number N of units, because 
the exponent l//x is greater than unity. Furthermore, Sn is not self- averaging but keeps 
on fluctuating in the large iV limit. The result (3.8) therefore appears as a fluctuating 
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scaling law. The expression of the distribution of the rescaled random variable £ M is 
obtained by using (3.9) in the first of the formulae (2.5), 

fcM= J %/ C '^- (3.10) 

Rey>0 

This probability density is referred to as the (properly normalized) Levy law of index /i. 
It is a non-trivial universal function in the sense that it is determined by one parameter 
only, namely the value of fi. References [16] provide a comprehensive presentation of 
Levy laws, and References [17] provide overviews of applications of Levy variables in 
Physics. 

Expanding e _yM as a power series in the integrand of (3.10) and integrating term 
by term with the help of the second of the identities (2.7), we obtain the convergent 
expansion 

Using the complement formula for the gamma function, we find the alternative form 

= £ r(i + M c-^\ (3.i2) 
k>i K - 7r 

The behaviour of the probability density /c (£) at large values of C is determined by 
the leading order term of this expansion, 

/^£)«^r(l + / ,)£-(^) (£-oo). (3.13) 

7T 

Its behaviour at small values of C is obtained from (3.10) by means of the method of 
steepest descent. We thus find an exponentially fast fall off: 

fcM -exp^l-ii)^)^ 1 "^ (£-0). (3.14) 

The density ferity admits a particularly simple expression in closed form for fj, — 1/2, 
which is 

„-l/(4£) 

fc y M = (3,5, 

The expression of the moment function rnc^ip) may be derived starting from (2.6). 
Using (3.9) in the right-hand side of the identity, we obtain it in three steps. We first 
evaluate the integral with the help of the relation ([14, p. 342]) 

i 

10 fi 

We then rewrite the result using the difference formula for the gamma function. We 
finally extend the obtained result to negative values of Rep by analytic continuation. 
We thus find the simple formula 

mcM - r r 1 (1 "_"^ > (Rej></<)- (3-17) 



f 

JO 
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3.3. Self-averaging with anomalous fluctuations: 1 < /i < 2 

In this case the distribution of the unit actions s n is still broad, but it has a finite 
mean (s). Using the second expression (2.11) of f s {x) in (3.1), we obtain 

f SN (x) « exp(iV (-(s)x + |r(l - fi)\ (soxY)) (3.18) 
in the relevant regime (N large and x small). Setting 

S N « iV(s) + s |r(l - ril^N 1 '" £ M , (3.19) 
we find that the moment generating function of the rescaled random variable £ M is now 

UM = ^- (3- 2 °) 

Equation (3.19) gives the scaling law of the sample action in the present case. This 
expression involves two different microscopic characteristics of the distribution f s (s), 
namely the mean (s) and the amplitude so of the power law tail (2.8). The fluctuations 
of S N about its mean N(s) increase as N 1 ^, with 1/2 < l/fi, < 1. On the one hand, 
these fluctuations are much smaller than the mean action in the large N limit, and 
so Sn is self-averaging. On the other hand, these fluctuations increase more rapidly 
than iV 1//2 , the power law of the central limit theorem involved in the expression (3.3). 
Furthermore, the distribution of the rescaled random variable £ M is not Gaussian. 

Using (3.20), we obtain that £ M is distributed according to the following Levy law 
of index ji: 

fc,(C)= J ^ vC+r - (3-21) 

Rey>0 

This probability density is universal, as it only depends on the value of the index /i. 
It is non-vanishing for all values of C Its fall off at large positive values of C can be 
derived by linearizing (3.21) as 

/r„0C)« / ^ e "V- ( 3 - 22 ) 

Rey>0 

Using again the second of the identities (2.7), we thus obtain 

fcM „ r(l + M ) (£^+oo). (3.23) 

7T 

This expression is the first term of an expansion analogous to (3.12). This expansion 
is, however, asymptotic but divergent in the present case. The behaviour of fc (C) at 
large negative values of £ is obtained from (3.21) by means of the method of steepest 
descent. We thus find 

-(/i-l)(Uj j (£^_oo). (3.24) 

Since the exponent /// (/i — 1) is larger than 2, this fall off is of a superexponential type. 
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4. Transmission statistics in the ensemble of samples with a fixed length L 

The statistics of the transmission probability through samples with a fixed length L are 
determined by the distribution of the sample action Sl, equation (2.2). The random 
character of Sl has two origins, as it is the sum of a random number Nl of independent 
random variables. We shall be interested here in the scaling properties of Sl in the limit 
of large values of L. 

It is again useful to introduce the moment generating function of the sample action. 
Its expression is now 

fs L (x,L) = '£p n (L)fs n (x), (4.1) 

n>0 

where p n {L) = Prob {Nl = n} is the probability that a sample of length L consists 
of exactly n units and fs n {x) the moment generating function corresponding to this 
particular number of units. Using (3.1), we obtain 

fs L {x,L) = Y.Pn{L)fs{x) n . (4.2) 

n>0 

When dealing with functions f(L) that depend on the sample length, it will prove useful 
to use their Laplace transform f(u) with respect to this variable, 

~ roc 

/(«)= / e~ uL f(L)dL. (4.3) 
Jo 

Equation (4.2) becomes 

W L (x,u) = J2Pn(u)f s (xT- (4-4) 

n>0 

The expression of the Laplace transform p n (u) can be obtained as follows. The 
definition (2.3) of Nl implies 

p n (L) = Prob {L n < L < L n+1 } 

= Prob {L > L n } - Prob {L > L n+1 } . (4.5) 
The Laplace transform of the first of these probabilities is simply 
^ Prob {L > L n } e~ uL dL = (j e - uL dL) 

= < e "" L ") = kM = (4. 6) 

u u u ' 

where the last equality has been obtained in analogy with the derivation of (3.1). The 
Laplace transform of the second probability is obtained upon replacing n by n + 1 
in (4.6). It follows that (see, e.g., [15]) 

p n (u) = Uu) n (n > 0). (4.7) 

Using (4.7) in (4.4), we find that the Laplace transform of the moment generating 
function has the following closed form expression: 

fs L M= l - f M (4.8) 

U(l - fl(u)f a (x)) 



Statistics of quantum transmission in one dimension with broad disorder 



11 



This equation is analogous to the Montroll- Weiss equation [18], which is well known in 
the theory of the continuous time random walks (see, e.g., [19]). 

In the following, the scaling law of Sl will be derived by using (4.8) in the small x 
and small u limits. The expressions (2.9) and (2.11) of f s (x) are different for fi > 1 and 
bi < 1, as well as those of fe(u) for v > 1 and v < 1. This leads to four different phases, 
which are identified by the conditions (/i > 1, v > 1), (/i < 1, v > 1), (/i > 1, v < 1), 
and (/j, < 1, v < 1). We shall discuss these four phases in this order. For the sake of 
simplicity we shall not distinguish the situations of normal and anomalous fluctuations 
when the sums are self-averaging. Said otherwise, we shall not discuss separately the 
cases fi > 2 and 1 < /i < 2, nor the cases v > 2 and 1 < v < 2. 

4-1. Phase I: fi> 1 and v > 1 

In this phase both means (s) and (£) are finite. The result (3.3) implies that S N iV(s) 
and similarly that Ln ~ N(£), up to fluctuations that become negligible in relative 
value in the large N limit. Eliminating N between both estimates, we obtain 

S l*J^L, (4.9) 

again up to fluctuations that become negligible in relative value at large L. In order to 
obtain more accurate estimates of the mean (Sl) and the fluctuations about it, more 
restrictive conditions on the distributions fg and f s are needed. 

Let us consider first the case of narrow distributions with /x > 2 and v > 2. Using 
the expressions (2.9) of f s (x) and fi{u) in (4.8) and expanding this formula to second 
order in x and u, we find that the mean and variance of the distribution of the sample 
action have the following expressions: 



((Sl)) » ^ULL^ULl. (4.11) 



(t) 2 ((s 2 ))±(s) 2 ((e 2 )) 

(if 

Equation (4.10) shows that the actual expression of the mean is not exactly the one 
given by (4.9). There is a finite correction term which is proportional to ((£ 2 )) — (£) 2 = 
(£ 2 ) — 2(£) 2 and so may be of either sign. Equation (4.11) indicates that the variance 
increases in direct proportion to L. The first correction to this leading order result, which 
has not been written explicitly, is a constant term for v > 3 but grows as (s) 2 L 3 ~ u for 
2 < v < 3. 

We also consider the case of narrow distributions with a complementary distribution 
function that falls off more rapidly than any power law. In this case Sl is self-averaging 
in the strong sense that each cumulant of its distribution increases in direct proportion 
to L: 

((S k L )) » a k L. (4.12) 
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The expressions of the coefficients a k can be obtained as follows. Equation (4.12) means 
that the moment generating function f$ L (x, L) has the following exponential growth for 
large values of the sample length: 

fs L (x,L)~e A ^ L . (4.13) 

Here A(x) is the generating function of the coefficients a k) 



A(x) = J2 L jra k . (4.14) 



k>i 



Equation (4.13) in turn implies that the Laplace transform fs L (x,u) has a simple pole 
of the form 

fa^-^wr (4 - 15) 

On the other hand, (4.8) shows that fs L (x, u) has poles whenever fe(u)f s (x) = 1. Hence 
the generating function A(x) is determined by the implicit equation 

fe(A(x))f s (x) = 1. (4.16) 

Using the expression (2.10) of f s (x) and fe{x) in this equation and expanding it as a 
power series in x, we can derive the expressions of the coefficients a k in a recursive way. 
We find 

a - {S) 

ai ~w 

a 2 = > ( 4 - 17 ) 

and so on. As expected, the expression of the coefficient a\ (resp. a 2 ) agrees with the 
one obtained from (4.10) (resp. (4.11)). 

4-2. Phase II: /i < 1 and v > 1 

In this phase the mean {£) is finite, but the mean (s) is infinite. The fact that (£) is 
finite implies that the number of units increases as N ^ L/{£), by analogy with (3.3). 
We therefore surmise that the scaling law of Sl is similar to the one of Sn in the case 
H < 1. 

The scaling law of Sl can be obtained as follows. Using the first expression (2.11) 
of f s (x) and the expression (2.9) of fe(u) in (4.8), we obtain 

(4.18) 

in the relevant regime (x and u small). The calculation of the inverse Laplace transform 
of fs L (x,u) gives 

fs L {x,L) » exp f-r(l - MsoxY^] . (4.19) 
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We note that this expression becomes identical to (3.7) if we replace the ratio L/(£) 
by the number N of units. This allows us to conclude that the scaling law of Sl in 
Phase II is 

S L ^s (T(l-^(^j ^ (4.20) 

where the rescaled random variable £ M is distributed according to the Levy law, 
equation (3.10). This result agrees with our surmise. 

4.3. Phase III: pL > 1 and v < 1 

In this phase the mean (s) is finite but the mean {£) is infinite. The fact that (s) is finite 
implies that Sn ~ N(s), by virtue of (3.3). We therefore surmise that the fluctuations 
of Sl are similar here to those of the number Nl of units in the samples. 

Following our surmise, we begin by deriving the scaling form of the expression of 
p n (L) in the limit of large values of L. Using the first expression (2.11) of fe(u) in (4.7), 
we obtain 

p n (u) « T(l - v)l^u v - x exp(-r(l - v)(e u) v n) (4.21) 
in the relevant regime (n large, u small), and so 

p n (L) « T(l - v)£l { p;u u - 1 exp(-T(l-is)(l oU yn + uL). (4.22) 

Reu>0 

Let us introduce a rescaled variable Z v by setting 



r 1-1/ 







This positive random variable has a limiting distribution which is easily obtained 
from (4.22). Setting z = uL, we find 

fzXZ)= I ^-V-^. (4.24) 



2vri 

Rez>0 

This probability density admits an expression in closed form for v — 1/2, case in which 
it is a half-Gaussian: 

e -z 2 /t 

/ Zl/a (Z) = — (Z>0). (4.25) 



"K 

The random variable Z v can be expressed as 

Z v = C v ~ v , (4.26) 

where L v is distributed according to the Levy law, equation (3.10). Indeed, setting 
Z = Cr v and z = xC in (4.24), we obtain the following expression for the probability 
density of C u : 

fcM = 2 J ' • (4-27) 

Rex>0 
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This expression is then put into the canonical form (3.10) by means of an integration 
by parts. 

The expression of the moment function of Z v is readily obtained by means of (3.17) 
and (4.26). We find 

mzM = WTW) V**p>-i). ( 4 - 28 ) 

This result implies that the distribution of Z v is narrow in the strong sense that all its 
moments are finite. Setting p = k in (4.28), we indeed obtain 

^ = WTW) < 4 ' 29 > 

for any integer k. 

The relations (4.23) and (4.26) may also be found in the following heuristic way. 
By analogy with (3.8), the sample length L N has the following fluctuating scaling law 
in the case v < 1: 

L N = £ (T(l-u)) 1 ^N 1 /' / C u . (4.30) 

Expressing N in this equation in terms of the other quantities, we obtain both above 
relations at once. 

Let us now turn to the scaling law of S L . Using the first expression (2.11) of fe(u) 
and the expression (2.9) of f s (x) in (4.8), we obtain 

» x £w%Lr (4 ' 31) 

in the relevant regime (x and u small), and so 

Is l (Sl,u) » r(1 ~ ( ^ X ' 1 exp (-r(l - vWouYj^j . (4.32) 

We note that this expression becomes identical to (4.21) if we replace the ratio S L /(s) 
by the number n of units. This allows us to conclude that the scaling law of Sl in 
Phase III is 

Our surmise has proved correct, as the same random variable Z u , equation (4.26), enters 



Z u . (4.33) 



the scaling laws of the number of units and of the sample action. 
4-4- Phase IV: /i < 1 and v < 1 

This last phase is the most interesting one, because neither the mean (s) nor the mean (£) 
is finite. This implies that Sn ~ N 1 ^ and Ljy ~ N 1 ^, by virtue of and analogy 
with (3.8) respectively. The combination of these two scaling laws gives S N ~ L v ^. We 
therefore surmise that the sample action follows a fluctuating scaling law of the form 
Sl ~ L"^ in this phase. 
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To obtain the scaling law of Sl, it is easier here to work with the Laplace transform 
ms L (p, u) of the moment function than with the one of the moment generating function. 
These two functions are linked together by the identity 

I poo 

ms L (-p,u) = — y ^ fs L (x,u)x p - 1 dx (Rep > 0), (4.34) 

which is readily obtained by taking the Laplace transform of (2.6) with respect to L. 
Using the first expression (2.11) of f s (x) and of fe(u) in (4.8), we find 

x r(i - 1/)^"- 1 , inrS 

fs ^ u) w r(i-,)(^r+r(i-^) (so ,r (435) 

in the relevant regime (x and u small). Putting this expression into (4.34), we obtain 
that of ms L (p, u) by following the same three-step procedure as in the derivation of the 
expression of the moment function m c ^(p), equation (3.17), the integral being evaluated 
with the help of the relation ([14, p. 292]) 

r°° , r(p/i/)r(i -p/u) , n „ x /A , 

/ dx = Ky/ p; v . Fl w (0<Rep<n). 4.36 

JO l + ZX* pLZ P ^ 

We thus find 

The inverse Laplace transform is straightforward, and leads to 

This result means that the scaling law of Sl in Phase IV is the following: 



r(i-v)J \e 

Here Y^ u is a positive rescaled random variable whose moment function is 

, . r(i +p//ji) r(i -p/n) / -p* \ ,„ N 
m ^ W= rii + Wrini-V (-/.<itop<M). (4-40) 

The comparison of this moment function with the one of the Levy variable £ M , 
equation (3.17), leads to the product relation 

m Y^ v (p) = m c^(p)m Cu (-pu/fj,) = m Cfl (p)m C[/ - v /^(p). (4.41) 

It follows that Y^ u can be expressed as 

V = £ M A,-^, (4.42) 

where the two independent random variables £ M and £„ are distributed according to 
Levy laws of respective indices \i and v. 

The expression of the scaling law, equation (4.39), and the relation (4.42) may also 
be found in a heuristic way, that is, by eliminating N between the fluctuating scaling 
laws (3.8) and (4.30). 
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The expression of the distribution of is obtained by using (4.40) in the second 
of the formulae (2.5), 

f Y (y)= 1^- m-p/rini+p/n) Y - P -i (4 43) 

JY » A > J 27ri T(l-p)T(l+pu/fi) { > 

This probability density is a universal function in the sense that it is determined by two 
parameters only, namely the values of the indices /i and v. 

The probability density fy (Y) can be expressed in the form of a convergent series 
of powers of Y^ in the case \x>v and inverse powers of Y^ in the case /i < v. In the 
former (resp. latter) case this series is obtained from (4.43) by summing the contributions 
of the poles of the integrand at p — —k/j, (resp. p = kfi), for k > 1. Using the difference 
and complement formulae for the gamma function, we find 

IrJY) = EM)" =^ 0. > »), 

k>i 71 L w , 4U) 

ky * T(i + ku) ^- ; 

The behaviour of the probability density at small (resp. large) values of Y is determined 
for all values of the indices \i and v by the leading order term of the above expansions: 

71 (4 45) 

Jy "" 1 ; 7T r(i + i/) v ; 

The power law fall off of fy v (Y) at large values of Y is similar to that of the Levy law 
at large values of C, as given by (3.13). Its prefactor hardly depends on u, as 1/T(1 +u) 
varies between 1, which is reached for v — or v — 1, and l/r(l + f m in) = 1.12917. . ., 
which is reached for the value z/ min = 0.46163 ... of z/ at which the function T(l + v) 
takes its minimum T(l + u m - m ) = 0.88560 . . . The behaviour of fy^iY) at small values 
of Y is different from that of the Levy law at small values of £, as given by (3.14). It 
is characterized by a power law divergence in Y~^~^\ whose amplitude is a decreasing 
function of v which vanishes in the v — > 1 limit. 

The preceding equations take a simpler form in the particular case of the symmetric 
situation where /i — v. Indeed, the scaling law (4.39) has then the form 

of a fluctuating scaling law growing proportionally to the sample length L. The 
relation (4.42) becomes 

Y m = % (4.47) 

Hence the random variable Y^ is distributed as the ratio of two independent random 
variables and C^, each one of them being distributed according to the same Levy law. 
Using the difference and complement formulae for the gamma function, (4.40) becomes 

™ Y Jp)= S ™ nP , ; (-// < Rep < //). (4.48) 



LY m (4.46) 
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The probability density of admits an expression in closed form. Indeed, each series 
in (4.44) can be summed and both give the same result, which is 

f Y (Y) = ^ (4 49) 

JY ^ { ] 7rY(Y» +y-^ + 2c0S7T/x)' 1 ' 

This probability density has been originally discovered by Lamperti [20]. 

The overall shape of the distribution of Y^ is better revealed by considering the 
random variable 

A^ = lnV. (4.50) 

Equation (4.45) implies that the probability density of A^ u falls off exponentially as 
A -> ±oo: 



, , . , sin7rz/ r(is) A 

/V (A) » =^ e^ A (A -> -oo , 

7T r(/i) 

sin7r/i T(l + /i) _„ A , 



(4.51) 



The moments of the distribution of A^ v may be calculated by taking advantage of the 
identity fh (p) = m Y^(—p)- Expanding the moment function of Y^,, equation (4.40), 
in powers of p, we find 

<A„„>=<^, «ApH (2 -^r V . (452) 

where C = 0.57721 ... is the Euler-Mascheroni constant, and so on. The first of these 
relations indicates that the mean of the distribution is positive (resp. negative) if /i < v 
(resp. /j, > u), the second that its variance decreases as the values of one or both indices 
increase. Mean and variance vanish as fi — > 1 and v — > 1 simultaneously, in agreement 
with the limiting result (4.59). 

If (j, — v, a closed form expression of the probability density of = InY^ is 
easily obtained from (4.49). We find the symmetric law 

f ( K \ — sin7r/i 

/A w (A)- 27r(coshM + cos7r/i) . (4.5dj 

The probability density /a m „(A) is represented in Figure 2 for ji — 0.6 and five 
different values of v. The curves result from a numerical integration of the following 
expression: 

1 f°° ye~ iyA 

/a„„(A) = — / — --p^ . w^r— , r dy. (4.54) 

2/x j-oo smh(7ry//i)r(l — iyjr(l + \yv j \i) 

This integral representation of the probability density is obtained by setting p = iy 
in (4.43) and using the relation T(l + kt)r(l — iu) = ttu/ sinh(vr-u). The figure confirms 
that the shape of /a m „(A) is very sensitive to the value of v on the negative side, but 
much less so on the positive side. The position of its maximum (shown by a heavy dot) 
moves only weakly as a function of u, and in a manner that is not monotonic. 

The scaling laws of the first three phases can be obtained from that of Phase IV by 
taking the right limits and using the proper identifications. 
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Figure 2. Graph of the probability density of the random variable A^ v = In 
for /j, = 0.6 and five different values of v. From top to bottom near A = 0: 
v = 1 (case in which Y^i = £ M is distributed according to the Levy law (3.10)), 
v = 0.8, v = 0.6 (thick line, case of the symmetric probability density (4.53)), 
v = 0.4, and v = 0.2. A heavy dot indicates the maximum of the probability 
density for each value of v. 



• For fj, — > 1, the results of Phase III are recovered with the identification 

(8) = Hm— . (4.55) 

In particular, — > 1, and so (4.42) becomes 

Yi„ = Z V = C-\ (4.56) 

in agreement with (4.26). 

• For v — > 1, the results of Phase II are recovered with the identification 

(£) = lim (4.57) 
In particular, £^ — > 1, and so (4.42) becomes 

= (4.58) 
in agreement with (4.20). 

• For /j, — > 1 and z/ — > 1 simultaneously, the results of Phase I are recovered with the 
identifications (4.55) and (4.57). In particular, (4.42) becomes 

Y n = 1. (4.59) 

The rescaled random variable becomes deterministic in this limiting case, which 
implies that Sl becomes self-averaging, in agreement with (4.9). 

Finally, we have addressed by means of numerical simulations the question of the 
convergence of the distribution of Sl for a finite sample length L towards the asymptotic 
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scaling law (4.39). We have generated ensembles of 10 8 samples for four values of L 
ranging from 10 2 to 10 5 . Each sample consisted of a sequence of pseudo-random values 
of £ n and s n , distributed according to the truncated power laws 

/.(*) = ^ (*>!)> m) = ^ ('>!)■ (4-60) 

These distributions have s — £ — 1. For definiteness, we have considered the symmetric 
case fx — v. In this case equation (4.39) becomes (4.46) which implies 

ln^A^. (4.61) 

The choice of the rather high value \x = v = 0.8 for the index was motivated by 
the expectation that finite-size effects are larger for larger /j, because of the presence 
of logarithmic corrections to scaling for fj, — 1. An example of such corrections is 
mentioned at the end of Section 5. Figure 3 shows histogram plots of our numerical 
data for ]h(Sl/L) for the four values of the sample length L. A slow and monotonic 
convergence of the data towards the asymptotic prediction (4.53) is observed. 




A 

Figure 3. Study of the convergence of the distribution of Sl towards the 
asymptotic scaling law (4.39). Symbols: histograms of numerical data for 
\u(Sl/L) for u = v = 0.8 and four values of the sample length L. Solid line: 
theoretical prediction for the asymptotic distribution of A w , equation (4.53). 



5. Discussion 

In the present work we have studied a model of quantum transport in one dimension 
based on a disordered array of units whose lengths £ n and actions s n are independent 
random variables with either narrow or broad distributions. We have successively 
investigated the statistical ensemble at fixed number of units (in Section 3) and the 
statistical ensemble at fixed sample length L (in Section 4). The model is richer in the 
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latter ensemble than in the former because the sample action Sl is the sum of a random 
rather than fixed number Nl of independent random variables. Equation (4.8) has been 
the starting point of the analysis in the ensemble at fixed sample length. This equation 
is analogous to the Montroll- Weiss equation [18], which is well-known in the theory 
of the continuous time random walks. Nevertheless the overlap between the present 
work and the abundant literature on continuous time random walks (see, e.g., [19]), 
and more generally on fractional diffusion (see, e.g., [21]), seems to be minute. This 
is essentially because our basic random variables (lengths and actions) are intrinsically 
positive whereas most works on fractional diffusion consider displacements which are 
symmetrically distributed random variables, and so all the results differ. Let us however 
mention that the distribution of our variable Y^„ appears in some form in Reference [22]. 

The main outcomes of this work, that is, those concerning the ensemble at fixed 
sample length derived in Section 4, are summarized in the phase diagram shown in 
Figure 4. Four different phases can be distinguished, depending on whether the indices \x 
and v are greater or less than unity. A fluctuating scaling behaviour for the action Sl, 
of the form 



has been shown to hold in Phases II, III, and IV. This fluctuating scaling law consists 
of three factors: a power law L a of the sample length, with a scaling exponent a whose 
value dictates the localization properties of the eigenstates, as detailed below; a rescaled 
random variable Y, which has a non-trivial asymptotic distribution with some universal 
character; and a non-universal prefactor a, which encompasses the microscopic details 
of the model. Table 1 gives the expressions of the exponent a and of the rescaled 
random variable Y in the four phases of the model. Remarkably enough, Y can always 
be expressed in terms of positive Levy variables £ M . 



S L ^aYL a , 




V 




L N «N<£> 




L N « N <l> 
S N «N<s> 
S L «L<s>/<£> 



1 



0, 




S N «N<s> 

L v 




v 







Figure 4. Phase diagram of the model in the plane of the indices u and v, 
showing the four phases labelled I to IV. The scaling behaviour of Sn and Ln 
at fixed N (Section 3) and of Sl at fixed L (Section 4) is given for each phase. 
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Phase 


Indices 


a 


Y 


I 


pL > 1, v > 1 


1 


1 


II 


fjL < 1, v > 1 






III 


// > 1, v < 1 


V 




IV 


//<l,i/<l 




v - r r ~ v lv 

1 [IV '"/J, t-v 



Table 1. Expressions of the exponent a and of the rescaled random variable Y involved 
in the fluctuating scaling law (5.1) for the sample action Sl in the four phases of the 
model. Phase I: Y = 1 is trivial, and so Sl is self-averaging. Phases II, III, and IV: Y 
is non-trivial, and so Sl is not self-averaging. 

The physical consequences of the fluctuating scaling law (5.1) are the following. As 
recalled in the introduction, in the usual one- dimensional localization problem, where the 
random potential has short-range correlations, the conductance, which is proportional 
to the transmission probability |T| 2 , is a strongly fluctuating quantity in the insulating 
regime of long enough samples (I > (). On the other hand, the action S is much better 
behaved, as it is self-averaging in a strong sense. Whenever the fluctuating scaling 
law (5.1) holds, the action itself is not self-averaging but keeps on fluctuating in samples 
whose length is arbitrarily large. The same feature holds for the sample-dependent 
effective localization length, defined as 

£l = ^, (5.2) 

for which (5.1) implies 

U » (5.3) 

The lack of self-averaging of the action Sl and of the effective localization length £ L , 
which occurs in three of the four phases of the model, is one of its most salient features. 
As recalled in the introduction, the lack of self-averaging of the action was also observed 
in the superlocalized regime of one-dimensional models with non-stationary random 
potentials whose fluctuations grow with distance [10, 11, 12]. We are therefore tempted 
to argue that a fluctuating scaling behaviour of the form (5.1) is a generic feature 
of quantum transport in a large universality class of one-dimensional systems with 
broad disorder, in the sense of random potentials with strong and long-range correlated 
fluctuations. 

Three different regimes of localization are allowed by the fluctuating scaling 
law (5.1), according to the value of the scaling exponent a. 

• a < 1 (Phase IV for /i > v and Phase III): Underlocalization. In this regime the 
typical sample action Sl grows less rapidly than linearly with the sample length L. 
The effective localization length therefore grows with L as the power law £l ~ L l ~ a . 
This underlocalized regime is intermediate between localized and extended. It is 
qualitatively reminiscent of the behaviour right at the critical point corresponding 
to the metal-insulator transition in higher dimension. 
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• a > 1 (Phase IV for /i < v and Phase II): Superlocalization. In this regime 
the typical sample action Sl grows more rapidly than linearly with the sample 
length L. The effective localization length therefore decreases with L as the power 
law £l ~ L"^ -1 ). The positive exponent a — 1 can be arbitrarily large if /i <C 1. 
This superlocalized regime is similar to that observed in one-dimensional models 
with non-stationary random potentials [10, 11, 12]. For such potentials, however, 
we have a — 1 + H/2, with the Hurst exponent of the random potential satisfying 
< H < 1, and so the value of a cannot be greater than 3/2. 

• a — 1 (Phase IV for /i — v): Fluctuating localization. In this regime the typical 
sample action Sl grows linearly with L, just as in the usual one-dimensional 
localization problem, where the random potential has short-range correlations. 
However, at variance with the usual situation, the action keeps on fluctuating in 
samples whose length is arbitrarily large. The effective localization length also 
keeps on fluctuating from sample to sample, without ever systematically increasing 
or decreasing as the sample length is increased. Its probability distribution is 
indeed asymptotically independent of the sample length. The term of fluctuating 
localization is therefore fully justified. Equation (4.46) leads to the expression 

& K A' (5 ' 4) 

where is distributed as the ratio of two Levy variables with the same index, 
equation (4.47). Therefore, apart from the prefactor £ /s , the effective localization 
length £l is also distributed as the ratio of two Levy variables with the same index /i, 
and so its probability density is also given by a Lamperti law (4.49). 

Along the borderlines which separate the various phases shown in Figure 4, that 
is, when at least one of the indices is equal to unity, the above results are affected by 
logarithmic corrections to scaling. The simplest case occurs in the ensemble at fixed 
number N of units. The approach of Section 3 leads to the following result in the 
marginal case /i — 1: 

S N « Ns (lnN + C + r]). (5.5) 

Here C is a non-universal constant which depends on the whole shape of the probability 
distribution f s {s) and r\ is a rescaled random variable whose universal probability 
distribution is characterized by the following moment generating function: 

Uv) = e ylnv - (5-6) 

The sample action Sjy may be called marginally self-averaging, as its mean grows as 
iVlniV whereas its fluctuations grow as N. 

Appendix. Scattering and transport in the weak transmission regime 



This appendix is devoted to the study of the scattering and transport properties of a 
single unit and of an array of N units. We first recall the transfer matrix formalism for a 
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single unit [3, 8, 9]. We then use this formalism to obtain the amplitude of transmission 
through an array of N units in the regime where the transmission through each unit is 
small. We derive the law of addition for the actions (A. 24), which is the starting point 
of the model investigated in this paper. 

A.l. Transfer matrix formalism for a single unit 

Let us first consider the scattering of a quantum particle by a single unit. The time- 
independent Schrodinger equation corresponding to this process is 

-ip"(x) + V(x)i/j(x) = Ei/j(x), (A.l) 

where the potential V(x) has an arbitrary form in the scattering unit (a < x < b) 
and vanishes elsewhere. In both half-lines (perfect leads) on either side of the unit, the 
wavefunction is a superposition of plane waves exp(±iA;x), with k = VE: 

( A e lkx + B e~ lkx on the left of the unit (x < a), (\2) 
\ C e ikx + D e~ ikx on the right of the unit (x > b). 

The amplitudes A, B, C, and D are linked together by linear relations of the form 

where M is the transfer matrix of the unit. Time reversal invariance implies that M 
can be parametrized as 

M=(; e {), (A.4) 

where the star denotes complex conjugation. In addition, conservation of the probability 
current leads to the condition that 

detM= |e| 2 - |/| 2 = 1. (A.5) 

Let r (resp. r') and t (resp. t') denote the reflection and transmission amplitudes 
for a particle coming from the left (resp. right). The first case corresponds to the values 
A — 1, B — r, C — t, and D — 0, the second to the values A — 0, B — t' , C — r', and 
D — 1. Thus we have 

1 f* f 

t = t=± r = - J -, r' = L (A.6) 

e e e 

These quantities satisfy the relations 

|r| 2 = |r'| 2 = 1 - |t| 2 , r't* + r*t = 0. (A.7) 

The above relations allow one to derive several equivalent parametrizations of the 
transfer matrix in terms of the reflection and transmission amplitudes; the one we 
choose here is the following: 
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A. 2. Many units in the weak transmission regime 

Let us now consider the transmission through a sample consisting of an arbitrary 
number N of units which are put end to end. The transfer matrix M. N of the whole 
sample is the ordered product of the transfer matrices M n of the units, 

M N = M n ---Ml (A.9) 

This equation is equivalent to the matrix recursion equation 

M N = M N M N -!. (A.10) 

By analogy with the expression (A. 8) of the transfer matrix for a single unit, we 
parametrize M. n as 

M-^i "X")- (A - n) 

where 1Z n and 7/v are the reflection and transmission amplitudes corresponding to the 
whole sample. Equation (A. 10) leads to the following non-linear recursion relations: 
J_ = 1 + ri/R^-i K N = r N + TZ N ^ 
T N t N T N ^\ t N T N _ ± T N tjv^/v-i tN^N-i 
with the initial conditions % = 1 and 7Zo = 0. Combining these two relations, we obtain 
a three-term non-linear recursion relation linking TJv to T^-i and T/v_2, which is 

Y~ = — + T^^T (e r ~T ¥~) ' < A13 » 

Setting 

r N = h^ls. (A. i4) 

JJ N 

and using the second relation of (A. 7), we find that is linked to £V_i and £V_2 by 
the following linear recursion relation: 

V N = V N ^ - rN ; lTN ' (V N ^ - |t w _ 1 | 2 P 7V _ 2 ) . (A.15) 

1 — \tN-l\ 

We are interested here in the expression of in the regime in which the probability 
of transmission through each unit is small. In this regime, neglecting |^tv~i 1 2 with respect 
to unity, (A.15) simplifies to 

V N ^(l-r' N _ irN )V N _ u (A.16) 

with T>i = 1. Iterating this equation, we obtain 

N-l 

V N ^ H(l-r' n r n+1 ). (A.17) 

n=l 

We are thus left with the following expression for the transmission amplitude 7/v in the 
weak transmission regime: 

N 

n t n 

T N ^ Trr ^ . (A.18) 

n=l 
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This leading order result is valid up to corrections which are proportional to the |t n | 2 
in relative value and could be derived from the full recursion relation (A. 15). The 
result (A. 18) can be recast in the following more appealing way: 

^11 1 , k V 

7at~£i- —h-, —t 3 ...t N -iz ; t N . (A.19) 

1 - r x r 2 1 - r 2 r 3 1 - r^r-jv 

Here each denominator corresponds to the geometric resummation of all the repeated 
reflections between two consecutive units. These are indeed the only internal reflections 
from the whole multiple scattering expansion that survive at leading order in the weak 
transmission regime. 

The leading order expression of the sample action is therefore 

N N-1 

S N = -\ In \T N \ 2 ~-|]Tm l^l 2 + Re E ln (! - < r n + i), (A.20) 

71=1 71=1 



N N-1 



i.e., 

Sn ~ E s n + Re ]T In (l - |r n ||r n+1 | e^^ 1 )) , (A.21) 

71=1 71=1 

where we have introduced the action s„ and the reflection angles 9 n and 6^ of each unit, 
according to 

s n = -\ In |t n | 2 , r n = \r n \e w \ r' n = \r n \e W ". (A.22) 

We now introduce our second hypothesis besides the weak transmission regime, 
namely that the internal reflection angles n and 9' n are random and uniformly 
distributed between and 27r. This hypothesis can be postulated as a mere simplifying 
assumption, in keeping with a long tradition in localization theory [4, 3, 6, 23]. It can 
alternatively be justified on physical grounds by considering that the reflection angles 
are rapidly varying functions of the incoming momentum k. This is especially true in 
the weak transmission regime, in which the length £ n of each unit is expected to be such 
that ki n 3> 1. As a consequence, any narrow distribution of incoming momentum will 
result in a uniform averaging over the reflection angles. Carrying out this averaging 
over the angles 6 n and 6' n in (A.21), we find that the second sum vanishes by virtue of 
the identity 



L 



2tt z k r2w 



\n(l-ze w )d9 = / e lfce d# = 0, (A.23) 

k>i k {2 

which holds for any number z such that \z\ < 1. 

We thus find that the action Sn of the whole sample is approximately given by 

N 

Sn~J2 s n- (A.24) 

71=1 

This simple formula, which can be referred to as the law of addition for the actions, is 
the starting point of the model investigated in this paper. 
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